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Abstract. During an adiabatic pumping cycle a conventional two barrier quantum device 
takes an electron from the left lead and ejects it to the right lead. Hence the pumped charge 
per cycle is naively expected to be Q < e. This zero order adiabatic point of view is in fact 
misleading. For a closed device we can get Q > e and even Q>e. In this paper a detailed 
analysis of the quantum pump operation is presented. Using the Kubo formula for the geometric 
conductance, and applying the Dirac chains picture, we derive practical estimates for Q. 



1. Introduction 

Understanding of charge transport in mesoscopic and molecular size devices is essential 
for the future realization of 'quantum circuits' [1]. Of particular interest are adiabatic 
processes that take electrons and move them one by one via a device. The simple 
minded peristaltic point of view of such process is misleading: such picture looks valid 
in the case of an open circuit, but breaks down once the pump is integrated into a 
closed circuit [2, 3, 4]. It is the purpose of this paper to further elaborate on the 
physics of quantum pumping in closed circuits, and to provide a detailed analysis of a 
prototype pump. The interest and the feasibility of realizing experiments with closed 
circuits is discussed in Section 1.3 of [5]. 

1.1. The prototype pumping device 

The prototype example for a quantum pumping device is the two barrier model (Fig.l). 
The one particle Hamiltonian is 

H{Xi{t),Xi(t)) = ^-p 2 + Xr(t)5(x - xi) + X 2 {t)S(x - x 2 ) (1) 
zm 

where m is the mass of the particle and (x,p) are the position and the momentum 
operators. The region x < x\ is the left lead, and the region x > x 2 is the right lead. 
We refer to the segment x\ < x < x 2 as the "dot region" . The driving is performed by 
changing X\ and X 2 in time. In an actual experiment the control parameters X\ and 
X 2 represent gate voltages. In order to talk about charge transfer we have to define 
also a current operator. In what follows we use the conventional definition 
g 

1 = — {pS(x - x ) + 5(x - x )p) (2) 
zm 

where e is the charge of the particle and xq is an arbitrary section point. The 
momentary current via different sections is in general not the same. But if we integrate 
it over a whole pumping cycle the result becomes independent of xq . 
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The pumping device can be used in two different configuration. In the case of an 
open geometry (Fig.l) the leads are attached to reservoirs that have the same chemical 
potentials. For simplicity one assumes a zero temperature Fermi occupation. 

In the case of a closed geometry (Fig.l) the leads are detached from the reservoirs 
and are connected together so as to have a ring. This means periodic boundary 
conditions over a large space interval (-L/2) < x < (L/2). Furthermore, the closed 
system is assumed to be strictly isolated from any environmental influences. The 
closed system can be regarded as a network with two nodes that are connected by two 
bonds one of length L D (dot region) and the other of length L w (wire region). The 
total length of the ring is L = L D + L w . We assume that L D <C L. 

1.2. Open (leads) geometry 

The open version of the two barrier model has been considered in Ref.[6] using the 
scattering matrix formalism of Biittikcr Pretrc and Thomas (BPT) [7, 8] . A typical 
pumping cycle is illustrated in Fig. 2c. In the 1st half of the cycle an electron is taken 
from the left lead into the dot region via the left barrier, while in the second half 
of the cycle an electron is transfered from the dot region to the right lead via the 
right barrier. Naively, by this peristaltic picture, it seems that at most one electron 
is pumped through the device per cycle. This expectation is supported by the formal 
calculation. Using the BPT formula one obtains [9] 

Q = (l-5 T )e (3) 

where < g T < 1 characterizes the transmission of the device during the charge 
transfer. In the limit g T — > 0, which is a pump with no leakage, indeed one gets 
Q = e. Otherwise one gets Q < e. 

1.3. Closed (ring) geometry 

Our interest is in the closed version of the two barrier model. A major observation is 
that the pumped charge Q is not "quantized" even if the device is closed and isolated 
from any environmental influences. Moreover, it can be larger than unit charge. In 
fact we can have Q ^> e. 

The analysis that we are going to present demonstrates and refines general results 
that were obtained in Rcf . [3] . There we have worked out an artificial circuit which has 
been modeled as a 3 site system. In the present paper we would like to work out a 
major prototype model that allows the desired comparison between results for closed 
circuit as opposed to that of open geometry [Eq.(3)]. 

We are going to use the Kubo approach to quantum pumping [2, 3, 4]. The "Dirac 
chains picture" which we further review in the next subsection makes a distinction 
between "near field" and "far field" pumping cycles. The near field result has been 
further considered in Ref.[10] using an extension of the BPT scattering approach to 
quantum pumping. 

1.4. The Dirac chains picture 

In order to analyze an adiabatic [11] pumping cycle we have first to understand the 
geometry of the parameter space. In fact the parameter space of the two barrier 
model is three dimensional (X\, X2, X3) where X3 = $ is the flux via the ring. In 
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practice we assume a planar <j> = pumping cycle, but for the theoretical discussion 
it is convenient to regard <I> as a free parameter. 

We ask what is the amount of charge which is transported via a section of the ring 
per cycle. For this purpose we have to calculate the current I = (X) at each moment. 
If we were changing the flux we would have by Ohm law / = — G 33< E> where G 33 is 
called the Ohmic (dissipative) conductance. But if we change (say) the parameter X\ 
then I = — G 31 X\, where G 31 is called the geometric (non-dissipative) conductance 
[12, 13, 14]. In general we can write dQ = -G 31 dXi - G 32 dX 2 and hence 

Q = jldt = <j>G-dX (4) 

where X = (X U X 2 ) and G = (G 31 , G 32 ). 

The elements of the conductance matrix G k ^ can be calculated using the Kubo 
formula. It turns out [11, 14] that in the adiabatic limit G 31 = B 2 and G 32 = —B\ 
where B is the "magnetic" field (2-form) which appears in the theory of Berry phase 
[11]. The sources of this field are Dirac monopoles that are located at the points 
of degeneracy. For the double barrier model a given level n can have a degeneracy 
provided X\ = X 2l and either $ = or <f> = irh/e modulo 2irh/e. In fact we have for 
each level (excluding the ground state) two Dirac chains of degeneracies as in Fig. 2d. 

From the above observation one easily draws the following conclusions: (i) We 
can get Q 3> e for a tight cycle around a Dirac chain if the degeneracy is in the 
pumping plane, (ii) We can get Q <C e for a tight cycle around a Dirac chain if the 
degeneracy is off the pumping plane. (iii) We can get Q <~ e for a cycle which is 
located in the far field of a Dirac chain. The existence of a far field region is not self 
evident. This constitutes a major motivation for the present study. 



1.5. outline 

The outline of this paper is as follows. In section 2 we clarify the starting point of 
the calculation, which is the Kubo formula. In section 3 we introduce a preliminary 
discussion of the expected results and their significance. Then in sections 4 to 11 
we analyze the pumping process in the two barrier model. In particular we find the 
dependence of Q on the "radius" of the pumping cycle, and make a distinction between 
"near field" and "far field" results. 



2. The Kubo formula 



Given a time dependent Hamiltonian Tt(X) with X = X(t) we define T = —dH/dX 
and would like to calculate the generalized conductance as defined by the expression 

(J) = = -G X (5) 

We label by n the adiabatic levels of the closed ring. The adiabatic states are defined 
by the equation 7i\n) = E n \n) with implicit X dependence. Using these notations the 
Kubo formula for the geometric conductance can be written as 

m{=j£n) 

The above expression assumes that only one level (n) is occupied. If several levels are 
occupied we have to sum over them. If we have more than one control variable, say 
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Xi and X 2 , then we have to use the more elaborated notation G 31 and G 32 in order 
to distinguish between different elements of the conductance matrix. 

If the pumping cycle crosses very close to a degeneracy, we can get from Eqs.(5-6) 
a very large current J, and upon integration we can find that the transported charge 
is Q S> e. In the next sections we shall develop actual estimates for Q. But first we 
would like to further illuminate the significance of Q. 

The Schrodinger equation can be written in the adiabatic basis, where the 
transformed Hamiltonian matrix takes the form H nm — E n S nm + XA nm , where E n 
are called the adiabatic energy levels, and A nm is a matrix that can be calculated 
using a well known recipe (which is summarized in section III of [3]). One regards 
X as the "small parameter" . If the system is prepared an instantaneous adiabatic 
state |n), then the instantaneous current is J = 0, whereas if it is prepared in an 
instantaneous steady state (an eigenstate of H nm ), then the instantaneous current is 
finite, say I = Iq. Accordingly one can question the physical relevance of Q: After all 
typically the initial preparation is an adiabatic state, and not a steady state. 

So let us consider an actual physical scenario. For simplicity we assume that 
only two adiabatic levels are involved: The occupied level n and next (empty) level 
m = n+ 1. To make the dynamical picture simple we use an analogy with the dynamics 
of a spin 1/2 particle, and consider the illustration in Fig. 3. We regard the state n as 
spin polarized in the z direction, and the state m as spin polarized in the — z direction. 
The instantaneous steady states of the Hamiltonian are polarized along an axis that 
has a small tilt relative to the z direction. 

Initially the spin is polarized in the z direction and therefore I(t = 0) = 0. For 
some time we have \X(t) — X(0)\ -C SX C where SX C is the relevant parametric scale 
for the variation of the adiabatic levels. During this time interval the tilt angle is 
approximately constant. The spin is performing a precession around the tilted axis. 
As a results we get I(t) ^ 0. In fact the maximum current is I(t) = 2J . We get this 
current after half period of precession. 

So we have a modulated current I(t) that equals upon averaging J . As long as 
\X(t) — X(0)\ -C SX C the precession goes on as described above. But on larger time 
scales we have to take into account the variation in the tilt angle. Consequently the 
modulation of the current is no longer in the range < I(t) < 2 Jo, but rather it is 
shifted and may increase. Still the average stays approximately Jo. 

Thus we see that in the actual physical scenario the average over I(t) is the 
same as the Jo of the instantaneous steady state. The validity conditions for this 
statement are essentially the validity conditions of linear response theory, which are 
further explained in Ref.[3]. The discussion above illuminates the justification for the 
use of the first order steady state solution of Kubo for the purpose of evaluating the 
pumped charge in an actual physical scenario. 

3. Charge transfer during an avoided crossing 

The pumped charge Q is obtained via the integral Eq.(4) using Eq.(6) for G. On 
the basis of the naive heuristic picture of the Introduction (and see Fig. 2) we expect 
that most of the contribution to the integral would come from a small segments in X 
space where the last occupied level has an avoided crossing with the first unoccupied 
level. Later we define precisely the region in X space where this assumption is a valid 
approximation. 
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Let us try to sketch what might come out from the calculation. Later on we are 
going to do a proper job. But before we dive into the detailed analysis (which is quite 
lengthy) it would be nice to gain some rough expectation. Given that our interest is 
focused in a small energy window such that E n ~ E wc define v E = {2E/m) 1 ' 2 . The 
mean level spacing in the energy range of interest is 

A = « E — (7) 

while the energy splitting A s at the avoided crossing might be much smaller. We 
define the following notations 

\Fnm\ = O-Q (8) 
\Inm\ = ~~j~\/9tp (9) 



A s /A = ^T~g~o (10) 

where both < g v < 1 and < go < 1 arc dimcnsionlcss. Note that go is related to 
the overall transmission ~g T of the device. The adiabaticity condition is 

A 2 

\X\ « —s- (11) 
and from Eq.(5) with (6) we get the current 



The time of the Landau-Zener transition via the avoided crossing is 

St « (A Q /a )/X (13) 
and hence the transported charge is 

Q « (1)5* « (l^) 7 e ( 14 ) 

where go and g v should be estimated in the region of the avoided crossing. We note 
that g v l (1 — go) is like the Thouless conductance, and can be regarded as a measure for 
the sensitivity of the energy levels to a test flux. We have pointed out and discussed 
this issue in Refs.[2, 3], and later it was derived [10] in the context of the scattering 
formalism. 

In Fig. 4 we display the numerically determined Q for various path segments. The 
horizontal axis is the \Xi — Xi\ distance of the path segment from the degeneracy 
point. As \X\ — X 2 \ becomes small go — > 1 and we get (3»e. But the asymptotic 
value Q m e which is observed for large \X\ — X 2 \ cannot be explained by such a 
simple minded calculation. A major objective of the detailed analysis is to illuminate 
the observed crossover. 



4. The model, basic equations 

The one particle Hamiltonian of the two barriers model depends on set of parameters 
(Ai,A 2 ,$). From now on we use units such that e = h = 1, and characterize the 
geometry by the dimensionless parameter 

b = L W /L D > 1 (15) 
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We write the wavefunction on the two bonds as: 

V> dot (z) = (2/L) 1/2 Vfe sin(^ D (x)) (16) 

V> wlro (x) = (2/L) 1/2 ^ sun> w (x)) (17) 

where (p(x) = kx + const. Given (X\, X 2 ) and assuming <f> = 0, the eigenstates of this 
Hamiltonian can be found by searching k n values for which the matching conditions 
on the log-derivatives are satisfied. This leads to the following system of equations: 

2m 

cot(v? W i) + cot(^ D1 ) = ~ x i (18) 
2m 

cot(<^ W2 ) + C0t(<£ D2 ) = ~ X 2 (19) 
<f D 2 - foi = kL D (20) 

<AV2 - <PW1 = (21) 

where it should be clear that <p D1 = ip D (xi) etc. The corresponding eigenenergies are 
E n — fc^/2m. A similar system of equation can be written for <I> = ir. We can find the 
Qw/Qd ratio for a given eigenstate via the matching condition on the wavefunction at 
either of the two nodes: 

^/q^ sin(iySw) = Vfc sin(<£ D ) (22) 

where (tp D ,ip w ) mean either (i^duVwi) °r (^2,^2)- The normalization condition 
implies that 

%fe + ^w«l (23) 

For an "ergodic state" we have q w 1 for both bonds. In general we characterize an 
eigenstate by a the mixing parameter 



such that 6 = means a definite dot state, while 6 — tt means a definite wire state. 
In practice we have < 9 < tt. 

In the numerical analysis we use units such that m = 1 and L D = 1. Given X\ 
and X2 we solve the system of equations for the k n and for the ip^ at the nodes. 
Then we determine q( n>> at each bond, and from the ratio we get 6^™^ as well. 

5. Regions in X space 

We focus on a small energy window such that the energy levels of interest are 
E < E n < E + dE. We characterize a barrier X8Q using its transmission 



g(x) 



1 + ( —x 



-1 



The last expression holds if g(X) <C 1. We can regard g\ = g{X\) and gi = g{X?) as 
an alternate way to specify X\ and X<i- The (Xi, X2, $) space is divided into various 
regions (Fig. 5a). There is a region where g\ and g 2 are of order one (|1 — g\ -C 1). 
There the delta functions at the nodes can be treated as a small perturbation. There 
is a region where 1/6 -C .91,52 <C 1. There each dot level mix with many wire levels. 
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This intermediate region will be discussed in a future work [15]. Finally there is the 
region in X space which is of interest in the present study: 

5i,32 « 1/6 (26) 

We shall argue that in this region the states are categorized into "wire states" and "dot 
state" , which mix only whenever the energy level of the dot "cross" an energy level of 
the wire (Fig.5bcd). This allows to use "two level" approximation in the analysis of 
the mixing. 

The degeneracies of the Hamiltonian occur at points (X {r \ iW) along the 
symmetry axis of X space. They are divided into two groups: those that are located 
in the <& = 0(mod(27r)) planes and those that are located in the $ = 7r(mod(27r)) 
planes. In Appendix A we find the explicit expression for X^ r \ It should be clear 
that each degeneracy point is duplicated mod(27r) in the $ direction, hence creating 
what we call a "Dirac chain" . 

There are only two non-trivial degeneracies in X space which are associated with 
a given level n. One is with the neighboring level "from above" and the other is with 
the neighboring level "from below" . 

Once we locate a degeneracy point we can make a distinction between "near 
field" and "far field" regions. The near field is defined as the region where we can 
use degenerate perturbation theory in order to figure out the splitting of the levels. 
In contrast to that, the far field is defined as the region where we can calculate the 
splitting of the levels by treating the dot- wire coupling as a first order perturbation. 

In Fig. 5a we show one pumping path that crosses in the near field region, and a 
second pumping path that crosses in the far field region. 



6. Eigenstate analysis 



In the theoretical analysis it is illuminating to map the behavior of k n and in 
(Xi, X 2 , $) space. See Figs. 6-7. It is not difficult to realize that the variation of O is 
bounded as follows: 



Vb 



1 

A 9 



+ 1/2 



< tan(0/2) < Vb 



1 

A 9 



-1/2 



(27) 



where g is either gi or g 2 . The derivation of this result is as follows: The matching 
conditions at a given node implies that (<£ D , f w ) are constraint to be on one of two 
branches which are illustrated in Fig. 8. With each point of a given branch we can 
associate a O value via Eq.(24) and either Eq.(18) or Eq.(19). It is a straightforward 
exercise to express O say as a function of (p D , and then to find its minimum and 
maximum values. Assuming that g <§C 1 one obtains Eq.(27). 

Given O^"- 1 we can extract what are the at each bond, and what are the ip^ 
at the nodes. By solving Eq.(23) with Eq.(24) we get 

q = 6(cos(e/2)) 2 (28) 
q w = (sin(0/2)) 2 (29) 

then from Eq.(22) with either Eq.(18) or Eq.(19) we find an expression for ip at a 
given node. In particular we get 

1/2 

(30) 



f 



1± -=tan(6/2) 
Vb 





'I 


) 


4 S . 
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Note that ip is well defined mod(7r). There are states which are dot-like (9 <C 1), and 
there are states which are wire- like (G ~ n). There are regions where a dot-like state 
mix with wire-like state leading to pair of states with 6^+' ~ 0( _ ) ~ n/2. From the 
above formula it follows that in the latter case 

1/2 

(31) 



^1 



(-) 



i 

I 5 



i 

7b 



i 

l 9 



1/2 



(tan(e (+) /2) + tan(e ( -V2)) 



(32) 



7. Eigenenergies and "mixing" from perturbation theory 

In order to find the splitting and the mixing of the two levels we use perturbation 
theory once in the far field analysis and once in the near field analysis. In both cases 
we use the following notations. The unperturbed basis is |i) with a dot- like state |1) 
and a wire-like |2). The perturbed eigenstates \n) are indicated as |+) and |— ). The 
Hamiltonian in both cases has the general form 

El ^ + ( Wu Wl2 
E 2 )^\W 2l W 22 

The eigenvectors are real so W\ 2 = W 2 i ■ The Hamiltonian can be written as a linear 
combination of Pauli matrices 

W 22 \ _ , (Ex -E2 W n - W 22 ' 



H = H + W = 



(33) 



H . E L +E 2 + Wn_ 



1 + 



we define 
A., 




W? 2 



tan(0) = 



2Wv 



(W n - W 22 ) + (Ex - E 2 ) 
The eigenenergies are 
'Ex 



E n — 



E 2 Wxx + W 22 



cr 3 + Wx2crx (34) 

(35) 
(36) 

(37) 



1+) 
I-) 



(38) 
(39) 



2 2 

The eigenstates are found by rotating a spin half around the y axis at an angle 9 
cos (6/2) 
sin (6/2) 

-sin (6>/2) 
cos (6/2) 

Assuming that the unperturbed basis consists of distinct dot-like and wire-like states, 
it follows that we can use the following approximation: 

9 (+) = 6 (40) 

e (_) = it -6 (41) 

Using Eq.(28) and Eq.(32) this implies 

b \ S m(9)\ (42) 



(+)„(-) 



2 
2 

7b 



i 

A 9 



1/2 



1 



|sin(0)| 



(43) 
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We note that in the strong mixing region we have 9 w tt/2. 

In the next sections we obtain explicit expressions for the "splitting" and the 
"mixing" using the above scheme. From the numerics (see e.g. Fig. 7) we see that 
these are in fact very satisfactory approximations. 



8. Far field analysis 



The eigenstates of the ring in the X region of interest as defined in Eq.(26) can be 
found using first order perturbation theory with respect to the wire-dot coupling. This 
is explained below. We have verified numerically that the approximation is remarkable 
unless we are very close to a degeneracy point. This defines our distinction between 
"far" and "near" field regions. In the latter case we present in the next section a 
complementary treatment using degenerate perturbation theory. 

The unperturbed Hamiltonian in the far field analysis corresponds to 
X\ = X 2 = oo. Using the notations as defined in the previous section we take |2) 
as the n th wire state, and |1) as the closest dot state from above. Hence 

2 



1 

2m 
1 

2m 

Using the formula 



Ei = 



Eo — 



(1 + [n/6] intogcr ) 



Wa = - 



1 



we get: 
\W l2 \ 



Am 2 X 

2L D 

vl 



[<9V> w PV (i) ] 



2L 



1 

~x[ 

1 



1 

X~ 2 
1 



X\ X2 



X\ X2 



(44) 



(45) 



(46) 



(47) 



(48) 



(49) 



2-\/L D L v 

where the ± sign in the expression for the dot-wire coupling depends on whether the 
dot and the wire states have the same parity or not. 

In case of a far field pumping cycle we start (say) with very high barriers, and 
then lower one of them, say X\. If we neglected the dot- wire coupling W12, the dot 
level would cross the wire level at a point X\ = X^ that can be determined from 
the equation E\ + Wn = E 2 - At the vicinity of the avoided crossing we obtain 



\W 12 \ 



L {bg (n)y/2 A 
Ztt 



where 



(50) 



(51) 



From the condition |Wi2| <C A we deduce Eq.(26) which defines our X region of 
interest. Furthermore, from the results of the previous section we obtain expressions 
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A, = bgW ±.\\ X -XW\\ (52) 

sin{9) = vimw^m (53) 

We use the following notation, which we regard as a measure for the distance of the 
pumping cycle from the nearby degeneracy: 



\X-XW\ 



t Xi - x "" f2+ l{i^) <54) 

The significance of this notation will be further clarified in the next section where we 
extend the analysis into the near field region. 

9. Near field analysis 

The unperturbed Hamiltonian in the near field analysis corresponds to 
X\ = X 2 = X^ r \ We can find a rough approximation for X^ using the analysis 
of the previous section, but in fact we can find the exact expression which is derived 
in Appendix A, where we also define the obvious notations k r and g( r \ For later 
calculation the following approximation is useful: 

X (r) „ _^£ = (55) 

V 9 ^ 

The unperturbed basis consists of the dot-like definite parity state |1), and the wire- 
like definite parity state |2). We recall that for these states G attains its extremal 
values as remarked at the end of Appendix A. The energies of the unperturbed states 
are 

E 1 = E 2 = E r = ±-h* (56) 
and the perturbation matrix is 

W n = bg^^(SX 1 +6X 2 ) (57) 

W 22 = g ( - r ^(5X 1 +5X 2 ) (58) 

\W 12 \ = Vbg^^SX^SX^ (59) 
where 

6X X = X 1 - X^ (60) 
SX 2 = I 2 -l' r > (61) 
Consequently we can determine both the splitting and the mixing of the levels: 

A s = bgW±-\\X-XW\\ (62) 
• ia\ ^ \Xi - X 2 \ 

sm{9) = Tb \\X-XM\\ ^ 
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In the above expression we have extended the interpretation of the distance measure 
as follows: 



\X - X^\\ = yj (X 1 +X 2 - 2X(r)) 2 + I\X 1 - X 2 



(64) 



Now we realize that the far field results of the previous section are formally a saturated 
version of the near field results with 



|*i - X, 



2|oo 



(65) 



10. The expressions for G 

We are now ready to calculate G from Eq.(6). There are of course G 31 and G 32 but the 
expressions look alike so we concentrate, say, on the case X = X\, and suppress the 
node indication subscript whenever possible. By definition T = —dTL/dX = 5(x — x\). 
The current operator has been defined in Eq.(2), but we still have the freedom to set 
xo as we want. So the natural choice for sake of calculation, is obviously x = x\. We 
shall expand later on the results that would be obtained if the current were measured 
via a different section. The matrix elements of the operators involved are 



2 
L 
e 
mL 



V 1 q^Q^ sin(^ n) ) sin(^ m) 

ft,;, 



/ (m) (n) \ k r , 

sin(vJf, -fn ) H 



. / (m) 



■ sin {ipl 



(66) 
) 



One should notice that T nm is real and symmetric with respect ton«m interchange, 
while T nm is antisymmetric and purely imaginary as implied by time reversal 
symmetry. Once we sum Eq.(6) over all the occupied levels, nm terms cancel with mn 
terms. Within the framework of the "two level approximation" the only remaining 
term involves the occupied level n and the next empty level m = n + 1 

ran I m [-^nm] 



G 31 (X 1 ,X 2 ) = 2- 



We recall the following expressions: 

b 



o (+) o (_) 



f D 



sm(6)\ 



lv4 +) -*4 _) l 



1 

2 

Vb 



1/2 



I 9 



1/2 



|sin(6>)| 



(67) 

(68) 
(69) 

(70) 



Upon substitution we realize that both in the near and in the far field we can neglect 
the second term in the expression for X nm . Consequently 

1 2ev E IXr-Xsl 



G 6L (X l7 X 2 ) 



1 / , n3/2 ev E 
l {9b) ~L? A 2 



sin(6»)| 



(71) 



b^g ||X-XM|| 3 

We observe that in the near field, where \\X — X( r )|| is essentially the Euclidean 
measure of distance, we get the field of a Dirac monopole as expected. But as we 
go to the far field the \Xi — X 2 \ contribution saturates as explained in the previous 
section. 
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11. The calculation of Q 

Having found G 31 and a similar expression for G 32 we can perform the integration of 
Eq.(4) in order to obtain Q. We already pointed out that most of the contribution 
in the regions of our interest comes from the avoided crossings. In the near field it is 
most convenient to make the integration along X\ — X 2 = iconst segments, while in 
case of the far field we make the integration along X 2 = 00 and X\ = 00 segments. 

It is not difficult to realize that in the near field calculation both segments 
(X 1 — X 2 = ±const) give the same contribution. This means that we simply have to 
do one of the integral and to multiply by 2. But in the far field one should be more 
careful. If we measure the current in the middle of the dot (as indeed done in the 
numerics) then the same rule applies. But if we measure the current (say) at node 1, 
then the predominant contribution to Q comes obviously from the X 2 = 00 segment, 
so the result should not be multiplied by 2. Performing a straightforward calculation 
of the dX integral, using 



and taking the above discussion into account, we get in the near field 



One can show that this result is in agreement with the rough estimate of Eq.(14). 
However, in the far field we have to substitute the saturated value of \X\ — X 2 \ leading 
to the result Q ~ e. 

The calculation in the far field does not care whether the pumping cycle encircle 
a $ = degeneracy or a <& = ir degeneracy. It is only in the near field of that we 
see the difference. This is clearly confirmed by the numerics (Fig. 4). But a closer look 
reveals that the far field numerical result for Q is somewhat smaller than 1. This might 
look like a contradiction with respect to the general expectations. The resolution of 
this puzzle is related to the limitation of the far field perturbative treatment. Within 
the framework of this treatment 0(™' changes from 0(") = n to 1 "' = as we lower 
(say) the Xi barrier. But in fact we know from section 8 that is bounded. This 
means that not all the the probability gets into the dot region. Consequently, if we 
integrate along the X 2 = 00 segment, we expect to get Q < e as observed. On the 
other hand, in case of a full pumping cycle, we have to cross from the X 2 = 00 segment 
to the Xi — 00 segment. This was neglected in our calculation. Thus if we had a full 
cycle we would expect to get Q w e in the far field as implied by the Dirac chains 
picture. 

12. Discussion 

We were able to derive an estimate for Q in the case where the pumping cycle is 
dominated by a single degeneracy. Within this framework we can still distinguish 
between near field (where Q » e) and far field (where Q ~ e) regions. Our results are 
in agreement with those of Ref.[3]. We note that an optional derivation of the near 
field limit has been introduced in Ref.[10] using an extension of the BPT scattering 
formalism. But the latter derivation was not suitable to reproduce the far field result 
because it has been assumed there that the charge cannot accumulate in the dot region. 




(72) 



Quantum pump in a closed circuit 



13 



It should be re-emphasized that the results that we have obtained assume that 
the pumping cycle is dominated by a single degeneracy. In a follow up work [15] we 
shall discuss the case where the charge transport involves many levels, such that the 
contribution of neighboring levels (in Eq.(6)) is negligible compared with contribution 
that comes from \m — n\ > 1 levels. 

The results that we obtain for a closed geometry are very different from those 
that are obtained for an open geometry. This is because the motion of the electron 
is "recycled" . In a more technical language this means that multiple rounds should 
be taken into account in the calculation of correlation functions. Refs.[16, 4] further 
discuss how the Kubo formula can be used in order to interpolate these two extreme 
circumstances. 

It should be clear that adiabatic transport becomes counter- intuitive if one adopts 
a misleading zero order point of view of the adiabatic process. Moreover, even within 
the "two level approximation" it would not be correct to say that Q is determined by 
peristaltic mechanism: namely, it is not correct to say that charge transfer is regulated 
by the Landau-Zener transitions. 

A peristaltic mechanism would imply Q ~ e. In the near field we have Q > e so 
we do not have such mechanism for sure. This is also reflected by having the same / 
at both nodes as discussed in the paragraph of Eq.(73). However, even in the far 
field, where the peristaltic picture seems natural, we have realized that it is an over 
simplification: also in the case of a far field cycle a finite fraction of Q is contributed 
during the intermediate stages of the pumping cycle. 

Appendix A. Finding the degeneracy points 

If system is symmetric (X\ = X2 = X) then we can distinguish between odd and even 
states leading to the following eigenvalue equations 



As we lower X we have an exact crossing whenever a dot state crosses a wire state 
with the opposite parity. We have an avoided crossing whenever a dot state tries to 
cross a wire state with the same parity. The later becomes an exact crossing if the 
flux through the ring is half integer. 

We can determine the degeneracy points by equating (A.l) with (A. 2). This gives 
sin(fcL w ) = — sin(fcL D ). For half integer flux it is convenient to use delta gauge on the 
middle of the wire, so as to get there ir phase jump boundary conditions. This implies 
that in the above equation we make the replacement (kL/2) 1— > (kL/2) + (n/2), hence 
getting the degeneracy condition sin(fcL w ) = +sin(fcL D ). We therefore conclude that 
we have degeneracies for 

k<r — ~ — rir (A-. 3) 

They are categorized into $=0 degeneracies for n r = 1, 3, 5, . . . and <!>=7r degeneracies 
for n r = 2, 4, 6, . . .. Their location is (X {r \ X (r ~>) where 



cot(fcL D /2) + cot(fc£ w /2) = X odd states 
tan(£;L D /2) + tan(fcL w /2) = X even states 



(A.1) 
(A.2) 



— cot(k r L D ) 



(A.4) 



m 



Accordingly 



g {r) = g{X (r) ) = sin 2 (£vL D ) = sin 2 (fc r L w ) 



(A.5) 
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At a degeneracy point the mixing parameter 9 that characterizes the odd and the 
even states attains the extremal values which are allowed by Eq.(27). This can be 
verified by deducing q m /q from Eq.(22) with ip v = k r L v /2 for the odd state and 
¥d = (tt/2) + k r L D /2 for the even state. 
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FIG. 1. Illustration of the model system. The two barrier pumping device is used in two different 
configurations. Left panel: open lead geometry; Right panel: closed ring geometry. The barriers are 
located at the nodes x\ and X2 while the current is measured through the section at xq. 




FIG 2. (a) Upper left: The energy levels of a ring with two barriers, at the beginning of the pumping 
cycle. It is assumed that the three lower levels are occupied, (b) Upper right: The adiabatic levels 
as a function of time during the pumping cycle, (c) Lower Left: The (X\, X'z) locations of the 
Dirac chains of the 3 occupied levels. Filled (hollow) circles imply that there is (no) monopole in 
the pumping plane. Note that for sake of illustration overlapping chains are displaced from each 
other. The pumping cycle encircles 2 + 1 Dirac chains that are associated with the 3rd and 2nd levels 
respectively, (d) Lower right: The 2 Dirac chains that are associated with the 3rd level. 
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FIG. 3. Cartoon of the adiabatic evolution 
within the framework of the two level 
approximation. The system is prepared 
in level n, and the nearby empty level is 
to = n + 1. The control parameter X(t) 
is being changed slowly, and therefore the 
system "oscillates" around the first order 
adiabatic solution. The energy of the latter 
is illustrated by a dashed line. Note that 
the identity of the adiabatic state changes 
gradually, and can be regarded as constant 



only on scales SX <C SX C . 




FIG. 4. Panels (a) and (b): Numerical calculation of the pumped charge Q. The 
model parameters are L u = m = e = 1 and L w = 3000.43. The current / is measured 
at the middle of the dot. The numerical integration is carried out along the segments 
which are indicated in panel (c), and the results are multiplied by 2 so as to include 
the equal contribution that comes from the second half of the cycle. There are two sets 
of data points. One set (filled circles) is for pumping cycles that encircle an in-plane 
degeneracy point (n r — 2993). A second set (hollow circles) is for pumping cycles 
that encircle an off-plane degeneracy point (n r = 2992). The location of the avoided 
crossing for each data set is indicated by the solid lines in panel (c). The near and 
that far field approximations that we derive for Q are indicated by the the solid lines 
in panel (a). The zoom in panel (b) reveals that Q in the far field is in fact slightly 
less then 1, which is explained in section 13. 
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FIG. 5. Regions in X space. The model parameters are the same as in the previous figure. Panel (a) 
displays the region of interest as defined in Eq. (26). It is bounded by the left and by the bottom solid 
lines which are defined by g(X) ~ 1/6. In-plane and off-plane degeneracy points are indicated by filled 
and hollow circles respectively. We indicate by arrow one in-plane degeneracy point (n r = 2993). A 
zoom of its near field is displayed in panel (b). The ellipse in panel (b) indicates a level splitting that 
equals A/10. The dashed lines in panels (a) and (b) indicate far field and near field pumping cycles 
respectively. In panels (c)-(e) we shows the energy levels (fcn) along three paths in X space, which 
are (Xi=X, X 2 =oo, $=0), and (Xi=X, X 2 =X, *=0), and (X\=X, X 2 =X, $=tt) respectively. The 
degeneracies n r = 2992. ..2999 are circled. The arrow indicates the representative degeneracy point 
n r = 2993. In panels (d) the odd states are indicated by dashed lines so as to distinguish them from 
the even states. 
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FIG. 6. The energy level splitting and the mixing parameter for two pairs of levels. The model 
parameters are L D = m = e = 1 and L w = 160.43. In panel (a) we show the contour lines for the 
energy level splitting of the first (even) dot level with an odd wire level (n = 158), and for the energy 
level splitting of the first (even) dot level with an even wire level (n = 157). The two cases are 
displayed again in panels (b) and (c) respectively where we plot both level splitting contours (solid 
lines) and contours (dashed lines). In the "even-odd crossing" case we have an in-plane degeneracy, 
which is indicated by a filled circle, while the inner most contour line is for A/5 splitting. Note that 
within the white regions the mixing is maximal (0 ~ it/2). In the "even-even avoided crossing" case 
the projection of the off-plane degeneracy point is indicated by a hollow circle. 
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FIG. 7. Tests of the perturbation theory based approximations (dashed lines) against the numerics 
(solid lines). The model parameters are L D = m = e = 1 and L w = 3000.43, and we focus on the 
degeneracy point n r = 2993. For these parameters X^ r ^ 465 and 4.5 X 10~ 5 . All the plots 

refer to the path (X2 — X\) = 5. In the left panel the dashed lines are derived from Eq.(37). In the 
middle panel the dashed lines are based on Eqs. (40-41) with 8 from Eq.(36). In the right panel the 
dashed lines are deduced from Eq. (30) 
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FIG. 8. The wire phase (^w/tt versus the dot phase ^ d /tt at a node with delta barrier 
g(X) = 0.225. The two branches are implied by the matching condition Eq.(18). The 
ratio |sin(</5 D )/sin(<p D )| attains its extremal values (Eq.(27)) at the points which are 
indicated by circles. 



